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Abstract 

This paper reports the initial results of a test series to evaluate a method for determining 
the normal incidence impedance of a locally reacting acoustically absorbing liner, located on 
the lower wall of a duct in a grazing incidence, multi-modal, nonprogressive acoustic wave 
environment without flow. This initial evaluation is accomplished by testing the methods’ 
ability to converge to the known normal incidence impedance of a solid steel plate, and to the 
normal incidence impedance of an absorbing test specimen whose impedance was measured 
in a conventional normal incidence tube. The method is shown to converge to the normal 
incident impedance values and thus to be an adequate tool for determining the impedance 
of specimens in a grazing incidence, multi-modal, nonprogressive acoustic wave environment 
for a broad range of source frequencies. 
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=height and length respectively, of the impedance tube test section 
=axial and transverse node numbers 

=number of points in the resistance and reactance grid, respectively 

=\Ta 

=free space wavenumber 

=the leading and trailing edge, respectively, of the test specimen 
=number of nodes in the x and y direction, respectively 
=two dimensional basis functions 
=number of upper wall measurement points 
=complex acoustic pressure at (x,y) 

=complex acoustic pressure at local nodes 1,2,3, and 4, respectively 
=reference pressure (20/iPa) 

=acoustic source and wall pressure, respectively 

=measured sound pressure in decibels 

= cartesian coordinates 

=locations of upper wall microphones 

=£ + ia, dimensionless admittance 

= global vector of the complex acoustic pressure 

= local vector of the complex acoustic pressure 

=measured phase angle for the complex acoustic pressure 

=9 A i\, dimensionless wall impedance 

= dimensionless exit impedance 

= local coordinate system for a finite element 


1 Introduction 

Highly efficient duct treatments for acoustic noise suppression continue to be a critical con- 
sideration in the achievement of environmentally acceptable commercial aircraft. To this end, 
a continuing concern in treatment technology is the accurate determination of the normal in- 
cidence impedance of an acoustic material underneath the boundary layers of grazing flows. 
Methods for determining the normal incidence impedance in such environments typically 
fall into three categories, “T-tube” method (refs. [1]), “in-situ” method (refs. [2, 3]), and 
the “propagation model” methods. All three methods have advantages and disadvantages, 
depending on facilities available, instrumentation and test liner construction or complexity. 
The disadvantages of the “T-tube” and “in-situ” methods are summarized in ref. [4], These 
two methods do, however, serve as useful complements to the “propagation model” method, 
which is the subject of this paper. 

The propagation model method is operationally convenient for extracting impedance 
spectra of test specimens in grazing flow environments. Consequently, flow duct facilities 
have been designed in an attempt to obtain the appropriate environment for duct propaga- 
tion models. The current procedure is to use an inhnite- wave-guide propagation model to 
extract the impedance of the test specimen from the measured data (i.e., spatial attenua- 
tion and phase rates on walls opposite the test liner) for a single, unidirectional propagating 
mode (refs. [4, 5, 6, 7]). However, in real facilities, these idealized test conditions (i.e., uni- 
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directional, single propagating mode) generally are attained only in an approximate manner 
under the best of laboratory conditions, and typically degrade severely with increasing mean 
flow. Additionally, many conventional test liner structures, and, in particular, the innova- 
tive liner structures currently under investigation, generate complex acoustic fields because 
of scattering of energy into higher order modes. While this may be a desirable result from the 
standpoint of achieving more efficient, broadband absorbing structures, it is a complicating 
feature that cannot be handled by the traditional unidirectional, single mode propagation 
model method. 

In a recent paper (ref. [8]), a method for extracting the acoustic impedance of a material 
installed as a finite length wall segment of a duct carrying a nonprogressive multi-modal 
sound held in the absence of how was presented. The method was shown to be successful 
for an assumed inhnite liner using analytically based input data. The purpose of this paper 
is to validate the impedance extraction method using measured data and a hnite length 
liner. More specifically, the impedance extraction method developed in ref. [8] is validated 
by testing the ability of the method to reproduce known normal incidence impedance data 
for a solid steel plate and a sound absorbing specimen. In contrast to the analytically based 
input data used in ref. [8], all input data used in this report were obtained from acoustic 
pressure measurements taken in a grazing incidence impedance tube facility with the test 
specimen installed. 

The remainder of this paper is organized into hve sections. Section 2 describes the 
measurement system used to obtain the input data for the numerical method and gives a 
brief description of each test specimen. The governing equations and boundary conditions are 
presented in section 3. Section 4 describes the numerical method which is used to extracts 
the unknown impedance. The known normal incidence impedances for the specimen was 
used as a baseline against which the accuracy of the model is judged. These results are 
presented in section 5. Finally, conclusions relevant to this paper as well as ongoing research 
activities are presented in section 6. 


2 Experimental Set-Up and Test Specimens 

The input data used to extract the impedance of each test specimen was obtained from 
measurements using a flow impedance tube in the NASA Langley Flow Impedance Labora- 
tory. This multi-configurational apparatus has a 51x51 mm cross-section and is designed to 
produce a controlled aeroacoustic environment with a Mach number of up to 0.6 over a test 
specimen length of 41 cm. Four 120-watt phase-matched acoustic drivers generate signals 
over a frequency range of 0.3 to 3.0 kHz, with sound pressure levels up to 155 dB at the test 
specimen leading edge. In the present investigation there was no mean flow. Sound pressure 
levels of at least 120 dB at each frequency of interest were set at the leading edge of the test 
specimen. A schematic of the flow impedance tube is provided in figure 1. The test section 
(the section of the duct between the source and exit plane) is 84 cm long. The upper and 
two side walls of the test section are stainless steel. 

Acoustic waves are propagated from left to right, across the surface of the test specimen, 
and into a termination section designed to minimize reflections over the frequency range of 
interest. Two 6 mm condenser-type microphones were flush-mounted in the test section, 
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one at a fixed location on the side wall, and the other on an axial traverse bar. A 13 mm 
wide precision-machined slot in the top wall of the flow impedance tube allows this axial 
traverse bar to traverse the test section length by means of a computer controlled digital 
stepping motor. The fixed location microphone was used to provide a reference for phase 
measurements. The data acquisition program automatically positioned the traversing mi- 
crophone at a number of pre-selected locations ranging from 18 cm upstream of the leading 
edge to 52 mm downstream of the trailing edge of the test specimen. At each measurement 
location a transfer function between the traversing and the fixed microphone was measured 
and used to determine the sound pressure level and phase relative to the fixed microphone. 
The complex acoustic pressure at the wall location is determined from the equation 

20 • i 20 

Pwall = Pref 10 SPL C * = p re f 10 SPL [ C OS </> + l SHI <p\ (1) 

A block diagram of the electronic instrumentation and signal conditioning system is given in 
figure 2. At each test frequency, a sound pressure level was set at the test specimen leading 
edge with a signal generator. 

The duct propagation model, which is discussed in the following section, requires a mea- 
surement of three sets of input data 

(1) The upper wall acoustic pressures 

( 2 ) The source plane acoustic pressure 

( 3 ) The exit plane impedance 

Unfortunately, the source plane acoustic pressure and exit plane impedance are functions 
of position along these planes. Therefore, transverse probe microphones should be used to 
measure this data when the test specimen is installed. However, a transverse probe micro- 
phone was not available for this initial investigation. The experiment was therefore carefully 
designed, so as to minimize higher order mode effects along the source and exit boundary. 
All data for the duct propagation model was therefore be obtained from measurements made 
by the upper wall traversing microphone. It should be noted that because of the sound ab- 
sorbing properties of the liner, it is not possible to avoid higher order mode effects in the 
liner region. In addition, higher order mode effects and reflections will occur in the vicinity 
of the leading and trailing edge of the specimen. 

In order to avoid the use of a transverse probe, the source plane was chosen 18 cm 
upstream (3.5 duct diameters) of the leading edge of the test specimen in the hardwall 
section of the duct, and the source frequency was kept below the cuton of higher order 
hardwall modes. Higher order mode effects caused by the installation of the test specimen 
are expected to decay upstream of the leading edge of the test specimen. Therefore, the 
source pressure at each point along the source plane was set to the value measured at the 
upper wall source location. 

A Similar procedure was applied at the exit plane. The exit plane was located 25 cm 
downstream (5 duct diameters) of the trailing edge of the test specimen, also in the hardwall 
section of the duct. The method developed in ref. [9] was used, with a hardwall test speci- 
men installed to obtain the exit impedance. This was done by measuring complex acoustic 
pressures at five locations along the top wall of the impedance tube using the axial traverse 
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bar mounted microphone. The same exit impedance was used for both the hardwall and the 
conventional liner configurations. Because the exit plane is 5 duct diameters downstream of 
the trailing edge of the liner and higher order modes are cutoff, higher order modes generated 
by the installation of the liner are not expected to carry appreciable acoustic energy to the 
exit plane. Thus, the exit impedance values at all points along the exit plane were set to 
that obtained at the upper wall. 

Input data were obtained for two specimens. The first specimen was a stainless steel 
plate (hardwall test specimen) for which the impedance was known. The second (a sound 
absorbing specimen) was a conventional uniform liner that consisted of a 100 MKS Rayl 
hbermetal facesheet bonded to a 100 mm deep aluminum hexcell honeycomb structure backed 
by a rigid plate. This type of liner construction has been used by the aeronautics industry for 
at least 35 years. Installed, this test specimen spanned the 51 mm width of the test section 
for a length of 41 cm. Normal incidence impedance measurements for this sound absorbing 
specimen were obtained using a conventional normal incidence impedance tube measurement, 
identical to that used in ref. [9]. Because of the simplification that the absence of mean flow 
affords, these measured normal incidence impedances and the impedances predicted using 
the method described in this report (using grazing incident sound) are expected to compare 
favorably. Thus, the measured normal incidence impedance is used as a baseline against 
which to judge the accuracy of the impedance extraction method when flow is absent. 


3 Governing Equation And Boundary Conditions 


In this section a brief description of the governing equation and boundary conditions used 
to extract the impedance of each test specimen is presented. Figure 3 depicts the applicable 
geometry and coordinate system used to model the test section of the Langley Impedance 
Tube. The axial and transverse directions are denoted by x and y, respectively. The spanwise 
direction, normal to the (x,y) plane, is not shown in the figure. To limit the propagation 
model to two dimensions, plane waves are assumed in the spanwise direction. This is a 
reasonable assumption since the sidewalls of the impedance tube are rigid and the sound 
frequency is below the cuton of any higher order modes. The lower and upper walls of the 
duct are located at y = 0 and y = H, respectively, and are also rigid except where the liner 
is located. The source and exit planes are located in the hardwall section of the duct at 
x = 0 and x = T, respectively. The test specimen (assumed locally reacting) is installed 
in the lower wall and has a normalized impedance, ( = 6 + fy. The leading and trailing 
edge of the specimen are at x = L\ and x = L 2} respectively. To perform the computation 
for the hardwall baseline (i.e. , the stainless steel plate) it is convenient to use the acoustic 
admittance of the specimen, defined as (3 = £ + io . Measurement points are located along 
the upper wall at x = aq, x 2 , x 3 . . ,x m as shown. Acoustic pressure measurements at these 
locations were obtained with the traversing microphone as described in the previous section. 

The mathematical problem is to find the solution to the Helmholtz equation 


d 2 p(x,y ) 
dx 2 


d 2 p(x,y ) 
dy 2 


+ k 2 p(x,y) 


0 


( 2 ) 


Along the source plane of the duct (x = 0), the sound pressure level and phase of the acoustic 
pressure held is measured and expressed as a complex pressure, p s , using equation (1). The 
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source plane boundary condition is therefore 


p(0,y) =p s 


( 3 ) 


The boundary conditions along the rigid upper wall and the rigid portion of the lower wall 
are equivalent to the requirement that the gradient of the acoustic pressure normal to the 
wall vanishes 


dp 

dy 


0 


( 4 ) 


At the exit plane (x 
the exit impedance 


L), the ratio of acoustic pressure to axial velocity is assumed to equal 


9p(L,y) 

dx 


~ ikp(L,y ) 

Cexit 


( 5 ) 


Throughout this work, all impedances are normalized with respect to the characteristic 
impedance, poC 0} of the air in the duct. Finally, the test specimen is assumed to be locally 
reacting, such that the boundary condition along the lower wall portion of the duct containing 
the test specimen is 

dp ikp 

ai = ~ (6) 


Equations (2)-(6) constitute a boundary value problem that can be solved to determine 
the unknown impedance, £, provided the upper wall pressure is known (ref. [8]). The solution 
for the unknown impedance can be put in terms of known functions only for special cases 
of the boundary conditions and upper wall pressure. These cases are not useful for general 
application. Therefore, the unknown impedance, £, must be determined by a numerical 
method for cases of practical interest. 


4 The Numerical Method 

The numerical method chosen to extract the unknown impedance of the lower wall is dis- 
cussed in ref. [8] and only sufficient detail is presented here for continuity and completeness. 
Because the experimental set-up was designed to minimize higher order mode effects, the 
foregoing discussion will assume a constant source pressure, p s , exit impedance, (exit, and 
lower wall impedance, (. However, these functions need not be constants to successfully 
apply the method (ref. [8]). When applied to the current acoustic problem, the method may 
be interpreted as an approximation of the continuous acoustic held as an assemblage of rect- 
angular finite elements as illustrated in figure 4. It is assumed that there are N nodes in the 
axial and M nodes in the transverse directions of the duct. A typical rectangular element, 
[/, J], is shown in figure 5. Each element consists of four local node numbers labeled 1, 2, 
3 and 4, respectively. Each element is considered to have width a = (xj +1 — xj) and height 
b = (yj+i —y.j) as shown. The objective of the method is to determine the impedance £, that 
minimizes the difference between the upper wall pressure obtained from the finite element 
solution, and that obtained from measured data with the specimen installed. 

To begin, Galerkin’s finite element method is employed to minimize the held error. The 
held error is distinct from the wall error function, which will be used to extract the unknown 
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impedance. First define the local coordinates, and rj for the element 


x 

a 

Now define the field error function as 



( 7 ) 


E(^ri) 


d 2 p(<^ } rj) 

ck 2 


+ 


d 2 p(<^ } rj) 
drf 


+ k 2 p(c ~, t]) 


Within each element p(<j, r/) is represented as linear combination of four functions, Afi 
and Afi which comprise a complete set of basis functions 


( 8 ) 

,N 2 ,N 3 


p(s,rj) = N 1 (q,rj)p 1 + N 2 (q,rj)p 2 + N 3 (q,rj)p 3 + N 4 (q,r])p 4 (9) 

Ni(^v) = fi(^)fi(v) N 2 (^r]) = f 2 ^)f 4 (r]) (10) 

N 3 (s,ri) = f 2 (q)f 2 (ri) N 4 (s,ri) = /i(?)/ 2 (??) (11) 

/i(^) = 1 — ^ = ^ (12) 


In an ideal sense, the solution to the sound field is obtained when the field error, E(<^ } r/) } 
is identically zero at each point of the domain. This is approximately achieved by requiring 
that the field error function be orthogonal to each basis function Ni(<^ } rj). Contributions to 
the minimization of the field error function from a typical element are 


ab [ [ E(q, rj)kdr] = (13) 

Jo Jo 

where is a 4x4 complex matrix for each element [/, «/], and a 4 X ]^ column 

vector containing the unknown acoustic pressure at each of the four nodes of the element. 
The coefficients in the local matrix, were computed in closed form. It should be 

noted that the second derivative terms in the field error function, E(<^ } r/) , were integrated 
by parts in order that linear basis functions could be used. Also, the wall and exit plane 
boundary conditions were substituted into these integrated terms. 

Assembly of the global equations for the computational domain is a basic procedure in 
the finite element method. Appropriate shifting of rows and columns is all that is required 
to add the local element matrix, directly into the global matrix, [A], Assembling the 

elements for the entire domain results in a matrix equation of the form 


m)m = m (i4) 

where [A(£)] is a complex matrix whose order is MN } and {$} and { F } are MNx 1 column 
vectors. The vector {$} contains the nodal values of the unknown acoustic pressure and 
{ F } is the zero vector until the pressure source condition is imposed. It is necessary to 
apply the source condition to equation (14) before a solution can be obtained. Satisfying 
the noise source boundary condition consists simply of setting all nodal values of acoustic 
source pressure at the source plane to the known value of source pressure, this introduces 
nonzero elements into the first M components of { F }. 

The global matrix [A(£)] generated by Galerkin’s Method following application of the 
source condition is a complex symmetric matrix. Fortunately, owing to the discretization 
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scheme used, it will also be block tridiagonal. The structure of the matrix [A(£)] after impos- 
ing source conditions is shown in figure 6, where the superscript T denote matrix transpose. 
Note that [A(£)] is a square block tridiagonal matrix whose order is MN. This global matrix 
contains a number of major blocks (A/, Bj) which are themselves square and block tridiago- 
nal as shown in the figure. The matrix elements in each major block are computed explicitly 
in terms of the yet unspecified impedance, (. Much practical importance arises from this 
block tridiagonal structure as it is convenient for minimizing storage and maximizing com- 
putational efficiency. Special matrix techniques exist for a solution of this structure 1 . All 
computation and storage is performed only on the elements within the bandwidth of the 
matrix [A(£)]. 

The unknown impedance, (, of the acoustic material is determined from the measured 
upper wall acoustic pressure. The procedure is to determine the impedance (, such that the 
pressure along the top wall reaches its’ measured value at each of the measurement points, 
xj. The procedure consists of repeatedly cycling through the solution to equation (14), 
obtaining a set of upper wall pressures for each impedance value. As each new set of wall 
pressures is computed, it is compared to the measured value until convergence is achieved. 

The idea is best illustrated as follows for a perforate. Define the unknown impedance as 

( = 9 + i X (15) 

where 6 is the resistance and x the reactance. Resistance values are positive whereas reac- 
tance values span the real axis 


0 < 6 < oo , — oo < y < oo (16) 

It should be apparent that searching the entire upper half plane of the resistance/reactance 
space for the unknown impedance is impractical. Therefore, introduce the tranformation 

y = cot(kd), 0 < kd < tt (17) 

and search for the unknown impedance in the (0, kd) plane, where 6 is limited to, 0 < 6 < 10, 
from practical considerations. 

The complex plane, (9 } kd) } is now divide into I MAX evenly spaced intervals in the 6 
direction and JMAX evenly spaced points in the kd direction, as shown in figure 7. The 
increment spacing A 9 and kAd are 

A9 = I MAX - 1 ’ kAd = JMAX - 1 
Thus, a point ( = ()/j in this uniform impedance grid is 

Cij = 0 i Ai X j , 9 I = (I-l)A9, X j = cot (J - l)kAd (19) 

If the measured upper wall pressures with the specimen installed are p wa ii{x n , H) } and those 
computed from the finite element solution with ( replaced by ()/j are p(x n} H) } then a measure 

1 Gaussian elimination with partial pivoting and equivalent row infinity norm scaling is used to reduce 
the rectangular system to upper triangular form. Back substitution is then employed to obtain the solution 
for the acoustic pressure at the MN node points 


( 18 ) 


8 



of the closeness of (ij to ( is given by the normalized wall error function, EW((u) 


1 m 

EW(Cu) = \ Pvjall{ x n, H) - p(x n , H) I (21) 

in which | | denotes the absolute value of a complex quantity, E max is the maximum value of 
EW for all points (u in the impedance grid, and m is the number of known wall pressures. 

Determining the unknown impedance of the material is now recast as a minimization 
problem. Thus, ( should be chosen such that EW(() is a global minimum. The global 
minimum is obtained using a two-step method. First, use a coarse grid in the impedance 
plane and tabulate the normalized wall error function to determine the location in that grid 
of an approximate minimum point (9 l ,kd l ). Next use a fine grid centered about (9 l ,kd l ) 
to obtain a new minimum. The location of the minimum point of this fine grid corresponds 
to the unknown impedance, ( of the material. Coarse grid calculations are performed with 
IMAX = 51, JMAX = 31, and A 9 = kAd = 0.1. Fine grid calculations are carried out 
using IMAX = J MAX = 21, and A 9 = kAd = 0.01. The coarse and fine grid discretization 
were determined from numerical experimentation in an earlier work (ref. [8]). 

5 Results 

A computer code implementing the impedance extraction method has been developed. The 
finite element matrix equation (14) is solved using a routine a highly developed software 
package and minimization of the normalized wall error function is performed internally by 
an in-house computer code. The unknown impedance, £, is returned by the in-house code. 
Results were computed using a Dec- Alpha work station and were not computationally in- 
tensive (i.e. , requiring only 0.5 seconds of CPU time for each point in the impedance grid). 
Analysis of the solid steel plate, for which the known admittance is zero (/3 = £ + eU = 0 + 0i) 
was conducted in the admittance plane for the sake of convenience. A 231x21 evenly spaced 
grid is used ( N = 231 and M = 21) in the finite element discretization for all calculations. 
This grid ensured that a minimum of ten elements per axial wavelength was used in the 
finite element discretization at the highest frequency of interest. 

Figure 8 shows a plot of the measured resistance and reactance of the normalized exit 
impedance as a function of frequency. Note that reflections are present in the test section 
since the resistance is not unity and the reactance has a non-zero value. Thus, the acoustic 
held in the Langley how impedance tube is not purely progressive, even for the hardwall 
test specimen. The measured resistance and reactance spectrum of the exit impedance is 
tabulated in the second and third columns respectively, of table 1. Resistance values range 
from a minimum of 0.87 at 3,000 Hz to a maximum of 1.13 at 2,600 Hz. Reactance values 
range from a minimum of -0.12 at2,800 Hz to a maximum of 0.12 at 2,300 Hz. 

Hardwall admittance predictions using the method proposed in this study are shown 
hrst. The predicted normalized conductances, £, for the hardwall test specimen, are shown 
in hgure 9. Predictions are shown for twenty-six frequencies ranging from 500 to 3,000 Hz, 



EW(Cu) = ^ (C/j) 

J-J m.n.cr. 
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in 100 Hz increments. The known and predicted conductance values are tabulated in the 
fourth and fifth columns, respectively, of table 1. The sixth column of table 1 is the error (the 
difference between the known and measured conductance values). The maximum error in 
the conductance prediction is only 0.02 (see tabulated results at 2,800 Hz). Figure 10 shows 
a graph of the predicted susceptance, a. The predicted and known susceptances, truncated 
to two decimal digits of accuracy, were zero at each frequency; therefore, tabular results are 
not shown. Overall, the predictions for the conductance and susceptance using the measured 
input data appear as accurate as those predicted in ref. [8] for analytically based input data, 
with a solid surface installed. 

An appropriate next step in the validation of the impedance extraction method is to 
install the conventional liner and predict the impedance of the soft test specimen. The liner 
and end effects will generate higher order modes and reflections in the test section that were 
not present for the solid steel plate. Figure 11 shows a comparison of the predicted and 
measured resistances, 0, for the conventional liner at six selected frequencies (500, 1,100, 
1,500, 2,000, 2,500 and 3,000 Hz). Predicted and measured resistance values for the conven- 
tional liner are tabulated in the second and third columns, respectively, of table 2. Predicted 
resistance values range from a minimum of 0.39 at 500 Hz to a maximum of 0.72 at 2,000 Hz. 
At each frequency, the predictions are close to the measured values. A maximum error of 
-0.10 is observed in the resistance prediction at 1,500 Hz (see the fourth column of table 2). 
Figure 12 shows predicted and measured reactance comparisons for the conventional liner. 
The data for this graph are tabulated in the fifth and sixth columns of table 2. The vari- 
ability in the reactance component of the impedance is noticeably greater than that of the 
resistance component. The maximum error in the predicted reactance is -0.20 and also oc- 
curs at 1,500 Hz (see column seven of table 2). Note that predicted reactance values are 
generally less accurate than the predicted resistance. 

Since the largest error in the predicted and measured normal incidence impedance for the 
soft specimen occurred at 1,500 Hz, each impedance value was used to predict the attenuation 
of the soft test specimen with the prediction program. The attenuation predicted with the 
measured normal incidence impedance (£ = ,63 + 2.48z) was 1.70dB and that obtained using 
the predicted impedance (£ = .53+2. 28z) was 1.67dB. Thus, for the chosen soft test specimen, 
the diffence between the predicted and measured normal incidence impedance (i.e. , .1 for 
resistance and -.2 for reactance) doesn’t affect the overall noise reduction peformance of the 
test specimen. 

Overall, resistance and reactance predictions with the liner installed match the measured 
normal incidence values quite well. However, the predictions are not as good as those ob- 
tained for the hardwall specimen. This is in contrast to the predictions obtained in ref. [8] 
with analytically based data, in which the predictions were just as accurate for soft as for 
rigid walls. There are several possible explanations as to why the predictions are less accurate 
with the liner installed when experimentally based input data is used. These are 

(1) Twenty elements across the duct may not be sufficient to resolve the transverse acoustic 

held in the neighborhood of the discontinuties at the leading and trailing edge of the 

test specimen. 

( 2 ) Installation of the liner generates some transverse dependence in the acoustic held at 

either the source or exit plane (these effects have been neglected in the predictions). 
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( 3 ) Differences between the predicted and measured values may be due to small errors in 
the measured normal incident impedance determination. 

A more likely explanation is that a combination of these effects is responsible for slightly 
less accurate predictions when the liner is installed. 


6 Conclusions and Discussion 

A duct propagation model for extracting the acoustic impedance of a test specimen using 
measured data has been tested for effectiveness. Although the method as presented here 
does not contain flow effects, it is extendable to mean flows with shear. Results of this 
study show that the method is extremely effective in extracting the impedance of hard and 
soft specimens for a range of source frequencies. There is only a slight degradation in the 
accuracy of the impedance prediction when high quality measured data is used instead of 
analytically based data. 

The results of this investigation represent a significant step in the development of a 
new and more complete impedance measurement technique. This method currently extends 
zero-flow impedance prediction capability to test specimens with non-uniform impedance dis- 
tributions in nonprogressive acoustic wavehelds, while maintaining impedance measurement 
accuracy. Additional validation experiments are planned, which should 

(1) validate the measurement capability for more complicated sources and at higher fre- 

quencies 

( 2 ) validate the applicability of the method to variable impedance liners 

Efforts are also underway to include the effects of shearing flows into the measurement 
technique. One important challenge in both the flow and zero-flow effort is to provide 
sufficiently accurate boundary condition data at the source and exit planes such as prescribed 
by equations (3) and (5). It will be recalled that in the source plane, the complex acoustic 
pressure distribution is required (and also the entrance impedance when flow is present) 
whereas in the exit plane the exit impedance distribution is required. 

Presumably, the successful measurement of the complex pressures at two closely spaced 
planes would permit acoustic particle velocities to be inferred, thus permitting calculation of 
the exit or entrance plane impedance. Although fraught with difficulty, this approach is being 
actively pursued by the development of in-flow/out-flow acoustic probes to directly measure 
acoustic pressure distributions in both the source and exit planes. A second approach to 
acquiring source and exit plane pressure distributions is also being pursued. In this approach, 
acoustic pressure measurements are confined entirely to walls of the source and exit sections 
of the flow duct. These pressures are then used to ’construct’ the inflow pressure and 
particle velocity fields from a modal decomposition of the hardwall pressure measurements. 
The advantages of this approach are 

(1) the avoidance of in-flow or out-flow measurements 

( 2 ) the intuitive appeal and physical insight provided by the modal propagation model 
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Frequency 
in Hz 

Measured 
exit impedance 

Conductance of the 
hardwall test specimen 

/ 

9 

X 

Known 

Predicted 

Error 

500 

1.03 

0.02 

0.00 

0.01 

0.01 

600 

0.97 

0.03 

0.00 

0.01 

0.01 

700 

1.05 

-0.01 

0.00 

0.01 

0.01 

800 

1.05 

-0.01 

0.00 

0.00 

0.00 

900 

1.03 

-0.06 

0.00 

0.00 

0.00 

1,000 

0.98 

-0.04 

0.00 

0.00 

0.00 

1,100 

0.98 

-0.03 

0.00 

0.00 

0.00 

1,200 

0.94 

0.02 

0.00 

0.00 

0.00 

1,300 

0.97 

0.07 

0.00 

0.00 

0.00 

1,400 

1.01 

0.05 

0.00 

0.00 

0.00 

1,500 

1.02 

0.06 

0.00 

0.00 

0.00 

1,600 

1.11 

0.01 

0.00 

0.01 

0.01 

1,700 

1.05 

-0.05 

0.00 

0.01 

0.01 

1,800 

1.04 

-0.06 

0.00 

0.00 

0.00 

1,900 

0.95 

-0.10 

0.00 

0.01 

0.01 

2,000 

0.92 

-0.05 

0.00 

0.00 

0.00 

2,100 

0.91 

0.02 

0.00 

0.00 

0.00 

2,200 

0.92 

0.06 

0.00 

0.00 

0.00 

2,300 

0.99 

0.12 

0.00 

0.00 

0.00 

2,400 

1.07 

0.07 

0.00 

0.00 

0.00 

2,500 

1.10 

0.08 

0.00 

0.01 

0.01 

2,600 

1.13 

-0.07 

0.00 

0.01 

0.01 

2,700 

1.05 

-0.11 

0.00 

0.01 

0.01 

2,800 

0.99 

-0.12 

0.00 

0.02 

0.02 

2,900 

0.90 

-0.10 

0.00 

0.01 

0.01 

3,000 

0.87 

0.00 

0.00 

0.01 

0.01 


Table 1 


Exit plane impedances and wall conductance values for the hardwall test specimen. 


Frequency 
in Hz 

Resistance of the 
conventional liner 

Reactance of the 
conventional liner 

/ 

Predicted 

Measured 

Error 

Predicted 

Measured 

Error 

500 

0.39 

0.44 

-0.05 

-0.55 

-0.66 

0.11 

1,100 

0.46 

0.40 

0.06 

-0.68 

-0.61 

-0.08 

1,500 

0.53 

0.63 

-0.10 

2.28 

2.48 

-0.20 

2,000 

0.72 

0.73 

-0.01 

-1.37 

-1.47 

0.10 

2,500 

0.53 

0.45 

0.08 

0.09 

0.17 

-0.08 

3,000 

0.42 

0.49 

-0.07 

1.35 

1.30 

0.05 


Table 2: Predicted and measured impedance for the conventional liner. 
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Figure 1: Schematic of the Langley Flow Impedance Tube. 



Figure 2: The electronic instrumentation and signal conditioning system for the Flow 
Impedance Tube. 
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ure 3: Modeled portion of Langley’s Flow Impedance Tube and coordinate system 
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Figure 4: Finite element discretization of two dimensional duct 
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Figure 5: A typical finite element, [/, J], and local node numbering system 
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Structure of the global stiffness matrix Structure of each major block, 

with major blocks, A/, Bi and Bj each x is a complex number 



Figure 6: Structure of the global matrix and major blocks 
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Figure 7: Impedance grid m the complex plane 
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Figure 8: Measured exit impedance for the Langley Flow Impedance Tube 
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Figure 10: Predicted susceptance for the hardwall test specimen 
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Figure 11: Comparison of measured and predicted resistance for the conventional liner. 
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Figure 12: Comparison of measured and predicted reactance for the conventional liner 
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